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Grover's algorithm for quantum searching is generalized to deal with arbitrary initial complex 
amplitude distributions. First order linear difference equations are found for the time evolution 
of the amplitudes of the marked and unmarked states. These equations are solved exactly. New 
expressions are derived for the optimal time of measurement and the maximal probability of success. 
| They are found to depend on the averages and variances of the initial amplitude distributions of 

■ the marked and unmarked states, but not on higher moments. Our results imply that Grover's 

^\ ' algorithm is robust against modest noise in the amplitude initialization procedure. 
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I '— ^ ' It is now firmly established that there exists a gap between the computational power of quantum and classical 
computers. A dramatic example of the speed-up offered by quantum computers is Grover's quantum search algorithm 
jl],^) for finding a marked element among N possible input values, in the presence of an oracle. On average a classical 
computer would need N/2 oracle- queries, whereas a quantum computer can accomplish the same task using merely 
CS| ' 0(y/N) queries. The importance of Grover's result stems from the fact that it proves the enhanced power of quantum 
computers compared to classical ones for a whole class of oracle-based problems, for which the bound on the efficiency 
t"*** 1 of classical algorithms is known. 

Grover's algorithm can be represented as searching a preimage of an oracle-computable boolean function, which 
can only be computed forward, but whose inverse cannot be directly computed. Such a function is F : D — > {0, 1} 
where D is a set of N domain values (or states) and the preimages of the value 1 are called the marked states. The 
problem is to identify one of the marked states, i.e., some v £ D such that F(v) = 1. Problems of this type are very 
0^ ' common. One important example, from cryptography, is searching for the key K of the Data Encryption Standard 
(DES) ||, gi ven a known plaintext P and its ciphertext C, where F = 1 if the pair of plaintext and ciphertext match 
[i.e., Ek{P) = C where Ek is the encryption function] and F = otherwise. Other examples are solutions of NP 
and NP-complete problems, which include virtually all the difficult computing problems in practice Q). 

A large number of results followed Grover's discovery. These results include a proof || that the algorithm is as 
efficient as theoretically possible || ; a variety of applications in which the algorithm is used in the solution of other 
problems |7Hf4| ; and recently, an experimental implementation using a nuclear magnetic resonance (NMR) quantum 
computer |15|. Several generalizations of Grover's original algorithm have been published, the first of which dealt 
with the case of more than one marked state |f6f| . The algorithm was further generalized by allowing an arbitrary 
(but constant) unitary transformation to take the place of the Hadamard transform in the original setting p7| . 

In this Rapid Communication, we generalize Grover's algorithm by allowing for an arbitrary complex initial ampli- 
tude distribution. We present an exact solution for the time evolution of the amplitudes under these general initial 
conditions. We find that the generalized search algorithm still requires 0(y/N/r) iterations, where r is the number 
of marked states, although the maximal success probability can be small for certain unfavorable initial amplitude 
distributions. The case of an arbitrary initial amplitude distribution is particularly relevant in the presence of unitary 
errors in the gates implementing the initialization step, such as over- or under-rotations. Such errors can result in a 
deviation from the uniform initial amplitude distribution assumed in the usual treatment of Grover's algorithm, and 
as detailed below, our analysis shows that the algorithm will still work in the presence of modest errors. 

We will now present the modified Grover algorithm and derive difference equations for the time evolution of the 
amplitudes in it. We then solve these equations exactly and analyze the results. Let k (t) [l(t)] denote the amplitude 
of the marked [unmarked] states after t iterations of the algorithm. It was shown in p6| that the amplitude of the 
marked states increases as: k(t) = s'm[uj(t + l/2)]/- v /r, where u> — 2 arcsin(-i/ r / N) . At the same time the amplitude 
of the unmarked states decreases as: l(t) = cos[ui(t + 1/2)]/%/ A — r. For N ^> r the optimal time to measure and 
complete the calculation is after T = 0(y/ N/r) iterations, when k(t) is maximal. In our modified algorithm we simply 
omit the initialization step from Grover's original algorithm. It thus consists of the following stages: 

I. Use any initial distribution of marked and unmarked states, e.g., the final state of any other quantum algorithm 
(do not initialize the system to the uniform distribution) . 
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2. Repeat the following steps T times: 

A. Rotate the marked states by a phase of ir radians. 

B. Rotate all states by ir radians around the average amplitude of all states. This is done by (i) Hadamard trans- 

forming every qubit; (ii) rotating the |0) state by a phase of 7r radians; (iii) again Hadamard transforming 
every qubit. 

3. Measure the resulting state. 

Next, we analyze the time evolution of the amplitudes in the modified algorithm with a total of N states. Let the 
marked amplitudes at time t be denoted by ki(t), i = 1, . . . ,r and the unmarked amplitudes by h(t), i = r + 1, . . . , N, 
where the initial distribution at t = is arbitrary. Without loss of generality we assume that the number of marked 
states satisfies 1 < r < N/2. Let the averages of the amplitudes be denoted by 



*(i) = i$>(*) 

for the marked states, and by 



r 



1 N 

i—r+l 

for the unmarked states. The key observation is that the entire dynamics dictated by Grover's algorithm can be 
described in full by the time-dependence of the averages. Let us define 

C(t) = ^[(N-r)l(t)-rk{t)]. (1) 

Consider any marked state fcj(t). In each step of the algorithm this state is flipped to fc^(i) = — fcj(i), so that the 
marked average becomes k'(t) = —k(t). The unmarked states, on the other hand, do not flip, so that the average 
over all states after the flip is: x(t) = jr[rk'(t) + (N — r) l(t)] — C(t)/2. Rotation by ir radians around the average 
is by definition: k[{t) -> 2x(t) - k'^t) and h(t) -> 2x(t) - k(t). Hence, h(t) -> C(i) + h(t) and k{t) -> C(t) - k(t). 
Therefore, the time evolution of all amplitudes (of both marked and unmarked states) is independent of the state 
index, and satisfies: 

ki{t+l) = C(t) + ki{t) i = l,...,r (2) 

k(t + l) = C(t)-k(t) i = r + l,...,N. (3) 

By averaging over the states in Eqs. (Q) and (Q) we find that the average marked and unmarked amplitudes obey first 
order linear coupled difference equations: 

k(t + l) = C(t) + k(t) (4) 
l(t + l) = C(t)-l(t). (5) 

These equations can be solved for k(t) and l(t), and along with the initial distribution this yields the exact solution for 
the dynamics of all amplitudes. We proceed to solve the recursion formulae for arbitrary complex initial conditions. 
Let: 



Using the recursion formulae (^) and (||) and a few steps of algebra employing the definition of C(f) given in ([!]), we 
find that f + (t + 1) = e iuJ f + (t) and + 1) = e" iw /_(i). Here u, which is real and satisfies 



cosw = f-2-^, (6) 
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is identical to the frequency found by Boyer et al. in fll6| . The time evolution can now be written as 

/+(*) = e*"7 + (0) 
/_(<) = e-*-*/-(0). 

Clearly, |/+(i)| and are time independent quantities. The average amplitudes are 



*(*) = -*y ^ [e iwt U(0) ~ e-^/-(0)] (7) 

r(<) = i[e fart / + (0)+e- iwt /-(0)]. (8) 

Together with Eqs. this provides the complete exact solution to the dynamics of the amplitudes in the 

generalized Grover algorithm, for arbitrary initial conditions. 

We turn to an analysis of several properties of the amplitudes and to a simplification of the result describing the 
dynamics. Let a and 4> (real or complex) be chosen such that a = y / /+(0)/_(0), and e 2 ^ = /+(0)//_(0). Using Eqs. 
(0) and (||), the average amplitudes can be expressed concisely as follows 



N-r 

k{t) = J asin(wf + 0) (9) 

f(t) = aco5(ut + <f>). (10) 

This shows that there is a ir/2 phase difference between the marked and unmarked amplitudes: when the average 
marked amplitude is maximal, the average unmarked amplitude is minimal, and vice versa. [Note that when the 
ratio F(0)/fc(0) is real, a and cj> become real, with a 2 = |Z(0)P + \k(0)\ 2 r/(N - r) and tan0 = \/r/(N - r)fc(0)/F(0)]. 
Subtracting Eq. (@) from E q- d), an d Eq. (|) from Eq. (|) one finds: 

ki(t + 1) - + 1) = fcj(t) - k(t) 
h(t + l)-l(t + l) = -[h(t)-l(t)]. 

This means that: 

Ah = fcj(0) - fc(0) (11) 
Ah = 1,(0) - F(0) (12) 

are constants of motion. This allows us to simplify the expression for the time dependence of the amplitudes: 

ki(t) = k(t) + Aki (13) 
h(t) = T{t) + (-If Ah, (14) 

where Aki and Ali are given by the initial amplitude distribution (at t = 0). 

Eqs. (^|), ([l0|), (|l3|) and (Q) are thus an alternative, simplified form describing the dynamics of the amplitudes. In 
this picture all marked states evolve in unison so it is sufficient to follow the time evolution of their average. The only 
feature distinguishing the states from one another is their initial deviation from the average. The same holds true for 
the unmarked states, up to an alternation about their average. 

From Eqs. (jl^) and (|l4|) it follows immediately that the variances 

<^(<) = ;!>(*) -*(*)l 2 ( 15 ) 
»=i 

i N 

vi{t) = jj— E ifc(*)-i(*)i 2 ( 16 ) 

i—r+l 

are time-independent. Now, when a measurement is performed at time t, the probability that a marked state will be 
obtained is P(t) = ££=i \h(t)\ 2 . Since all the operators used are unitary, the amplitudes satisfy the normalization 
condition: 
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E^(*)i 2 + E i^(«)i 2 = i 

2—1 ?__r+l 

at all times. Using (y — y) 2 = y 2 — y 2 {y is a random variable), we find from Eqs. ( |l5| ) and (p"6|): 

r 

J2Mt)\ 2 =ral+r\k(t)f 

i=i 

AT 

E i^wi 2 = (^-^ 2 + (^-^)i^)i 2 - 

z=r+l 

Therefore, the probability of measuring a marked state at time t is given by 

P(t) = P av - APcos2[u;* + Re(0)] (17) 

where 

P „ = 1 - (JV - r)of - \ [{N - r) \l(0)\ 2 + r |fc(0)| 21 



(N - r)l(0) 2 + rk(0) 2 



The maximal value that this probability can obtain during the evolution of the algorithm is 

P — P 4- AP 

1 max — J av ~ • 

Given an arbitrary initial distributions of r marked and N — r unmarked states, with known averages k(0) and 1(0) 
respectively, the optimal measurement times are after 

T = [(j + 1/2)tt - Re(0)]/w 

iterations, for j = 0,1,2,... when the probability of obtaining a marked state is P max - An important conclusion is 
that to determine the optimal measurement times, all one needs to know are the average initial amplitudes and the 
number of marked states. Expanding u> in the expression for T in r/N <C 1 (at j = 0) one finds that the number of 
iterations before the optimal measurement probability P max is obtained is 0(y/N/r). However, the value of P max can 
vary significantly, depending on the statistical properties (average and variance) of the initial amplitude distribution. 
The expected number of repetitions of the entire algorithm until a marked state is obtained is 1/P max . 

We next consider the shapes generated in the complex plane during the time evolution of k(t) — {Re[fc(f)], Im[fc(t)]} 
and 1(f) = {Re[Z(i)], Im[Z(i)]}. Eqs. (H) and (|io| ) turn out to be identical to the equations that describe the polarization 
of electromagnetic plane waves [M . By this analogy, the contours generated by these equations are ellipses in the 
complex plane. The major axis of the ellipse of l(t) subtends an angle r\ with the real axis, where rj is given by 
e"' = The length of the major [minor] axis of the ellipse is a = \a\ cosh(Im</>) [b = \a\ sinh(lm0)]. Here a and 

4> are the parameters which appear in Eqs. (^) and jic|). The ellipse of k(t) has a similar shape, but its major axis 
subtends an angle 77 + 7r/2 with the real axis and its major and minor axes are longer by a factor of v/ (N — r) jr. 

When the ratio l(Q)/k(0) is real, one can easily show that |/ + (0)| = |/_(0)|. In this case the amplitudes evolve along 
a straight line in the complex plane, in analogy to the case of linear polarization of light, and P max = 1 — (N — r)af. 
The best case, in which P max = 1, is obtained for Grover's original (uniform amplitudes) initialization, where of = 0. 

A limit in which the algorithm is totally useless is obtained when either /+(0) = or /-(0) = 0. In this case 
the success probability P(t) remains constant during the evolution of the algorithm. This corresponds to the case of 
circular polarization of light. The worst case appears when P max = P(t) — for any t. In this case, which is obtained 
when <7? = /+(0) = /-(0) = k(0) — 1(0) — and (N — r)af = 1, the algorithm can never find the marked states. 

Finally, consider the case where the average and variance of the initial amplitude distribution are not known, but 
different runs of the algorithm use initial amplitudes drawn from the same distribution. Naively, one could pick a 
random number of iterations T r and thus find a marked state with probability P(T r ). Correspondingly, the expected 
number of repetitions of the entire algorithm using the same T r would be 1/ P(T r ) until a marked state is found. 
However, P(T r ) could be very small. A better strategy is now shown. From Eqs. (^|) and ([l?]) it follows that the 
period of oscillation of P(t) depends only on r/N, while the details of the initial amplitude distribution are all in the 
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phase 4>. Consider the case where one runs the algorithm twice, taking measurements at times T\ and T2 respectively, 
where T2 — T\ = w/(2w). From Eq. ( |l7|) it is clear that in one of the two measurements P(T) > P av > P max /2. In 
this case we need twice as many repetitions to obtain at least half the success probability compared to the case when 
the optimal measurement time is known. The slowdown is thus at most a factor of four. 

In this work we generalized Grover's quantum search algorithm to apply for initial input distributions which are 
non-uniform. In fact, it was shown that by simply omitting the first step of Grover's original algorithm, wherein 
a uniform superposition is created over all elements, a more general algorithm results which applies to arbitrary 
initial distributions. To analyze the algorithm, we found that the time evolution of the amplitudes of the marked 
and unmarked states can be described by first-order linear difference equations with some special properties. The 
most important of these is that all amplitudes essentially evolve uniformly, with the dynamics being determined 
completely by the average amplitudes. This observation allowed us to find an exact solution for the time-evolution of 
the amplitudes. A significant conclusion from this solution is that generically the generalized algorithm also has an 
0(^J N/r) running time, thus being more powerful than any classical algorithm designed to solve the same task. An 
important future application of these results is in the study of the robustness of Grover's algorithm against errors in 
the unitary operations used to implement the algorithm. Our results imply that the algorithm can tolerate a moderate 
amount of noise in the amplitude initialization procedure. Work extending these results to the case of errors in the 
inversion about average step, and in the case of an arbitrary unitary transformation, is in progress. 

This work was initiated during the Elsag-Bailey - I.S.I. Foundation research meeting on quantum computation in 
1997, and was completed during the ensuing meeting in 1998. 
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